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Abstract 

We study 1 + 1-dimensional SU(N) gauge theories with adjoint scalar matter representations, 
based on a dimensional truncation of 2 + 1 and 3 + 1-dimensional pure QCD, which approximate 
the dynamics of transversely polarized gluons. The glueballs are investigated non-perturbatively 
using light-front quantisation, detailed spectra and wavefunctions being obtained for the large- N 
limit. In general there is some qualitative agreement of the spectra with lattice Monte Carlo data 
from the higher dimensional QCD. From the light-front wavefunctions we calculate (polarized) 
structure functions and interpret the gluon and spin content of glueballs. We discuss the phase 
structure of the reduced theories in relation to matrix models for relativistic non-critical strings. 
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1 Introduction 



While there is little doubt that Quantum Chromodynamics |T[ provides a correct description of 
most high energy processes inside hadrons, our ability to determine its predictions for low energy 
properties, such as spectrum and wavefunctions, is however much weaker. Of critical interest are 
the hadrons with excited gluonic modes, glueballs, and hybrid mesons, whose existence is currently 
amenable to concerted experimental investigation for perhaps the first time. QCD has often been 
studied in 1+1 dimensions, especially in the generalization to an infinite number of colours N — * oo, 
beginning with 't Hooft The reasons for doing this are both kinematical and dynamical. Phase 
space in two dimensions is smaller and therefore simpler. Also many of the characteristic features 
of hadrons, which are known or expected to arise from highly quantum dynamical effects in 3 + 1- 
dimensional QCD (QCD 4 ), materialize as classical or mildly quantum properties in QCD 2 . Since 
there is no angular momentum in one space dimension, conspicuously absent, however, are any 
features of QCD 4 which rely on transverse degrees of freedom. In particular there are no gluons at 
all. Is it therefore possible to embellish QCD 2 so as to include some (discrete) remnant of transverse 
fields without unduly complicating the theory? In this paper we will study in light-front formalism 
the glueball boundstates of a 1 + 1-dimensional gauge theory based upon a (classical) dimensional 
reduction of pure QCD 4 , a model of gluons which retains all polarizations of these particles but 
restricts their 3-momenta to the form k = (k, 0,0). The result is QCD 2 coupled to two mutually 
interacting adjoint scalars, corresponding to the two possible physical polarizations of (massless) 
gluons in four dimensions. Despite the rather crude approximation of degrees of freedom, our 
analysis in general can yield qualitative agreement with spectra obtained by lattice Monte Carlo 
simulations of pure QCD 4 . Having said this, our aim here is not to reproduce the numbers of 
higher dimensional gauge theory exactly, since the reduced theory is not well-suited to this task, 
but to treat the longitudinal and some of the transverse dynamics precisely. In this respect the 
non-abelian reduced gauge model in 1 + 1-dimensions seems to capture much of the physics of the 
higher dimensional theory. 

In the late seventies a remarkable attempt was made by Bardeen et. al || to compute the glue- 
ball spectrum directly from a strong coupling expansion of large- iV light-front Hamiltonian lattice 
gauge theory.]] However, the expansion scheme ultimately left an effective potential for transverse 
motion which had to be determined empirically. Our calculation, based upon the dimensional 
reduction mentioned above, is formally similar, but obviously does not include transverse kine- 

1 Review of more recent approaches to solving four-dimensional light-front gauge theories with further references can 
be found for example in ref.B. 



2 



matics, which would require knowledge of this effective potential. In order to partially overcome 
the effects of this omission it is expedient to allow certainly one, and possibly two, adjustable 
parameters in the reduced theory, which are allowed by 1 + 1-dimensional gauge invariance as we 
shall see. The analysis we perform is a natural extension of more recent work on the spectrum 
of adjoint matter in 1 + 1-dimensional large- iV gauge theories M. The simpler dimensionally re- 
duced model of QCD 3 at large- N has been studied in ref.||, questions relevant to the Hagedorn 
behaviour (of adjoint fermions) in refs.Q, and the heavy particle potential in ref.||. In this paper 
we will first perform in more detail the broad spectrum analysis of refs.|5|, |(| in order to elucidate 
the implications for glueballs in QCD 3 , for which there is some preliminary lattice data ||. 

Our computations are performed in the large-iV limit, which still however allows the number 
of gluons to fluctuate, unlike the case for quarks. This is done both for practical reasons, in 
order to keep the Hilbert space small enough to be manipulated on a workstation when studied 
numerically, and also possibly for reasons of principle in order to avoid subtleties of vacuum 
structure. The approach we use is analogous to the original one of 't Hooft [|j i.e. naive light-front 
quantisation in light-front gauge. While this is not in general good enough for studying vacuum 
properties, such as the value of condensates, it appears to be sufficient for the excitation spectrum.^] 
To solve the boundstate problem numerically we employ discretised light-front quantisation [ fill 
12] conveniently implemented in terms of a Lanczos algorithm, described in an appendix. We 



check the numerical results analytically in a phase space wavefunction (PSW) approximation. 
The solutions yield a mass spectrum and also light-front wavefunctions from which virtually any 
structure function of interest can be extracted. 

The organisation of the paper is as follows. Section 2 introduces the light-front dimensionally 
reduced model. New numerical and analytic solutions of the boundstate equations for the reduced 
model of QCD 3 are described in section 3 and compared with the available lattice data. Following 
the physical lessons learned from this, a similar analysis is applied to the reduced model of QCD 4 
in section 4; the glueball boundstate spectrum is compared with lattice data and the structure 
functions are investigated. In section 5 we briefly look at the phase structure of the reduced model 
in more detail, particularly in relation to matrix models of string theory which have been studied 
in a similar light-front formalism [O, 14]. Finally, section 6 summarises the physical conclusions 



we draw from this work about non-perturbative QCD. Descriptions of the numerical procedures 
are as far as possible collected in the appendix. An extension of the results to mesons and hybrids 
is treated in [001. 



2 For a careful treatment of these questions see for example M, 
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2 Dimensional Reduction. 



We start from SU(N) Yang- Mills theory in D dimensions 



1 



F, 



4 g2 J T ^ F " U 

d p A v — d v A p + i[j4^, Ay\ , 



(1) 
(2) 



and choose light-front co-ordinates x^ 
this case 

5y = — ;r [ d D ~ 2 xdx + dx~Tr 
g 2 J 



(x° ± x D 1 )/y / 2 and light-front gauge A_ = 0. In 



-a + A p d-AP + \{d-A + f + A + J + - \F prj F^ 
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i[A p ,a^A%-d^d p A^; i,j = l,...,N ; p, a = 1, . . . , D - 2 . 



(3) 
(4) 



Repeated indices are summed. The field A + is a constrained variable which does not propagate in 
light-front time x + , leaving only the transverse gluons A p as physical degrees of freedom. Much 
of the complexity of gauge theory is due to the linear term in the longitudinal momentum current 
(HJ), so one might as an approximation study the theory truncated to zero modes d p A p = only. 
This is the same as beginning with a 1 + 1-dimensional adjoint gauge theory 

2'" 4g 2 



Sr 



dx + dx Tr 



; a, (3 = ± 



(5) 



where g 2 = g 2 / J d D ~ 2 x, D a = d a + i[A a , .], (j) p = A p /g, and t = —1. Sr inherits a subgroup of the 
Poincare and gauge symmetries of the D-dimensional theory. The residual internal symmetries of 
Sr resulting from gauge invariance of S are just the 1 + 1-dimensional gauge symmetries 



A a -» UA a U ] + \{d a U)U ] , 



Ucf> p U^ 



(6) 



Due to the asymmetry of the dimensional reduction procedure, which singles out one spatial 
direction, it is appropriate to allow different longitudinal and transverse gauge couplings. Hence 
we take t in (J^) to be a free parameter, which is still consistent with the reduced gauge invariance 
(|6|). In fact the potential for transverse fields (ft p in Sr is the only quartic gauge invariant term 
with the centre symmetry <j) p — > <j) p + const. . Higher order interactions are not suppressed 
on dimensional grounds, but neither are they required for renormalisability; in the interests of 
simplicity we will therefore use the action Sr. It acquires the 50(1, 1) subgroup of SO(D — 1, 1) 
Lorentz symmetries of S and also has discrete symmetries: charge conjugation C {A^- — > —A^) 
induces 

C : 4>ij -> -(fiji (7) 
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on the 1 + 1-scalars; there also remains the one-dimensional parity 

Pi : x D - x -> -x D ~ l ; (8) 

finally there are the internal symmetries under exchanges of the p index which, together with Pi, 
form a remnant of the SO(D — 1) spatial rotations. 

Analysing the boundstates of Sr by light-front quantisation, and taking the large- N limit with 
g 2 N fixed, results in additional simplifications. The calculation is similar to, but much richer than, 
the one performed by 't Hooft for the quark-anti-quark pair Q]. The calculation is richer since, 
unlike the fundamental representation for quarks, any number of adjoint particles may be involved 
in the formation of a boundstate. With x + as canonical time, the A + field can be eliminated by 
its constraint equation to yield canonical momenta: 



P+ = - 

t 

2 



J dx-Ti[d-<]) p d-(j) p ] ; (9) 



P~ = — I dx~Tr 



1 » j+±r 



(10) 



2 vrpi Ycriiv > r j <-> Tp 

and the zero-charge condition / dx~ J + = 0, where J + = — i[<f>p, d-<j) p \. As in the t'Hooft model, 
the latter condition restricts one to the space of singlet states under ^"-independent gauge trans- 
formations, whose masses are gauge invariant. Quantising on the null line x + = 0, 

[(j>pij{x~),d^ ak i{x~)) = ^i5 pa 5(x~ -x~)5ii5 jk , (11) 

the boundstate problem reduces to the solution of the mass shell condition, 

2P + p-\^ >= M 2 \^ > , (12) 

which, for a suitable basis diagonal in P + , is a Schrodinger equation for the light-front hamiltonian 
P~. Such a basis is formed by the fourier oscillator modes aj,; 

0Mx + =0)=*r J*L (Wfc+)e- ifc+ *" + al 3l (k + y k+x -) , (13) 



2vr Jo ^2k 

Kv(k + ),al lk Ck + )] = 5 pa 5(k + - ~k+)5 u S jk (14) 

These create the light-front Fock space by action on a particle-free vacuum a p (fc + )|0 >= 0. We 
will often drop the superscript on k + for clarity and refer to it as momentum (it equals k° in the 
traditional infinite momentum frame). Due to the zero charge constraint and the large- N limit 
the light-front wavefunctions one needs to consider are the set of singlet states of the form 

oo r p+ 



l*0P + )>=E I dk l ...dk n 5\ [ Y / km-P + j l p 2 --' " J Tr[at i (fc 1 )---at ti (fc w )]|0> . 

(15) 



n=2" u \m=l 
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For SU(N), Tr[0] = 0, and the one-gluon state is absent. The glueballs (15) may be pictured as a 



superposition of rings of flux created by a closed chain of n gluons with momenta fei, . . . , k n . The 
states with more than one Trace, which are decoupled from these stable flux rings when 1/N — > 0, 
correspond to multiple glueball combinations and may be neglected. Explicit oscillator expressions 
may easily be obtained for ||. 

The main objective of this paper is to solve (12) for D = 3 and D = 4, and to explore the 



consequences. The elementary interaction processes occuring at order g 2 are shown in Figure 1, 
the principal difference between D = 3 and 4 being the absence of the four-gluon transverse- 
magnetic i-terms in the former case (fig. 1(c)). We shall see that this difference has significant 
consequences for the accuracy of valence gluon approximations and the phase structure of the 
theory. In the large-iV limit these proceseses must form Feynman diagrams that can be drawn 
on a cylinder, whose ends are the flux rings (|l5|), rather than a strip as in the case of quarks. 
The only renormalisation which needs to be dealt with in the two-dimensional theory concerns the 
self-energies shown in fig. 1(a). There is no transverse gauge invariance to prevent the generation 
of masses for the fields 4> p and so a mass counter-term must be added to Sr if masslessness were 
to be maintained at the quantum level. In fact, as in the t'Hooft model, the most divergent part 
of the self-energy is cancelled by zero momentum A + -exchange (E-term in fig. 1(b)) when two or 
more gluons are present. We elaborate these observations in the specific cases to follow. 

3 Reduced QCD 3 . 

3.1 Boundstate Equations. 

First we reconsider the reduced model for D = 3 ||, ||. There is only one transverse gluon field, so 
we drop the p index and call f pl — Pn = f n - By substituting the oscillator expressions ([[D into P ± , 
normal ordering and taking matrix elements, the mass shell condition (|l^) can be conveniently 
summarised by an infinite set of coupled Bethe-Salpeter integral equations for the coefficients 
f n (xi, ■ ■ ■ x n ) of n-gluon Fock states (cyclically symmetric in their arguments); 

tTt/ti = 9 A7 fn{xi,X2,---iX n )-\-- fn{ x l, X2, ■ ■ ■ , X n ) 

g z f\ g z J>J x\ 4yxiX2 

+ / —, . {E[x 1 ,X2,y}fn(x 1 ,x 2 ,,...,x n ) 

J0 ^X\X2V\X\ + X2 - y) 

+ (D[x 1 ,x 2 ,y] - E[x 1 ,x 2 ,y})fn(y,x 1 +x 2 -y,.. -,x n )} 
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+ ay dz f n+2 (y,z,xi - y - z,x 2 , . . . ,x n ) 

jo jo <±^Jx\yz{x\ — y — z) 

F[x 3 ,x 2 ,x 1 ] + F[x 1 ,x 2 ,x 3 ] 

+ - , j == r f n - 2 {Xi+X 2 + X 3 ,X4,...,X n ) 

4Vxix 2 x 3 (xi + x 2 + xs) 

+ cyclic permutations of (xi, x 2 , . . . , x n ) , (16) 
where the functions 

PI i (xi + y)(xi + 2x 2 -y) 

£?[xi,x 2 ,y] = (gi _ y)2 (17) 

D[xx,x 2 ,y] = E[xi,-y,-xi] (18) 

F[x 3 ,x 2 ,xi] = E[x 3 ,xi,-x 2 ] (19) 



correspond to the elementary processes of fig.l. In (16) we have switched to the Lorentz-invariant 
momentum fractions x m = k^ n /P + , which are the Bjorken scaling variables; the two-dimensional 
light-front gauge theory exhibits exact scaling. The particle 4> has logarithmic and linearly diver- 
gent self-energy in 1 + 1 dimensions as a result of the first diagram in fig. 1(a). This diagram arises 
from normal ordering the J.J term in P~ ([!(]) (the currents J themselves are assumed already 
normal ordered) and gives an induced mass Q 

1 2 g 2 N [°° dp g 2 N f k k 

As in the t'Hooft model, the linearly divergent part of m 2 nd cancels against zero-momentum A + - 
exchange (fig. 1(b)) between two particles and can be absorbed into the .E'-term in ([!(]) to make 
the integrals manifestly finite at the Coulomb pole x\ = y. To the rest of m 2 nd we have added a 
bare </>-mass in order to leave a finite renormalised mass m in (|l6|). In this paper we will mostly be 
interested in the massless limit m — >• 0, although no pathology would be incurred in choosing any 
m 2 > 0. In this way one can be sure that whatever constituent gluon masses appear to manifest 
themselves in the spectrum are dynamically generated. (We stress however that agreement with 
the higher dimensional lattice results can be improved if a constituent mass m is put in by hand.) 
Note that dimensional reduction of classical pure QCD 3 can be achieved by compactification of 
one direction to very small radius, but at the quantum level this results in a diverging Debye mass 
for Polyakov loops [p]]. At small radius these loops can be identified with ^-fields and therefore 
our mass renormalisation prescription is entirely different from the highly compact gauge theory, 
whose spectrum is quite disimilar to the uncompactified one [O] . 

The elementary processes (fig. 1) for adjoint particles represented in ( JT6[ ) include the Coulomb 
potential (E integral), which is linear in two dimensions; it would also be present for fundamental 
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representation particles but with half the strength g 2 — > g 2 /2 since an adjoint source has two flux 
lines attached to it rather than one. In addition in ( |l6| ) there is a 2 — > 2 annihilation channel (I? 
integral) and pair creation and annihilation of gluons (F integral), which are not suppressed by 
the large-iV limit but for low-lying levels are found to be dynamically suppressed to a very high 
degree ||. 

A survey of the broad features of the mass spectrum M 2 resulting from large- N adjoint scalars 
in two dimensions has been performed in refs. |3|, ||. A massive spectrum of stable glueballs is 
organised into approximate valence gluon trajectories for low-lying levels, while at higher ener- 
gies pair creation dominates to give a complicated picture. We have concentrated on the highly 
structured seven or so lowest states resulting from (16) in order to compare with Monte Carlo 
results of the (unreduced) lattice QCD 3 . Simple approximate analytic solutions may be found to 
the equations (|0]) while numerical solutions may be straightforwardly obtained with the help of 
Mathematica and a workstation. By cutting the interval of allowed momentum fractions < x < 1 
into x £ (1/K,3/K,5/K, ■ ■ ■) for some integer K, the problem becomes one of finite matrix di- 
agonalisation (Discretised Light-Cone Quantisation) [17|| , the continuum limit being achieved by 
extrapolating K — > oo. We found it useful to employ a Lanczos algorithm in computing and 
diagonalising the mass matrix M 2 in this way. A Hilbert space of dimension ~ O(10 3 ) can be 
comfortably handled. Details of our numerical work are given in the appendix. In addition to mea- 
suring the masses M from extrapolation of finite K calculations to the continuum limit K = oo, 
we also used the (normalised) light-front wavefunctions |^ > to compute structure functions of 
these mass eigenstates at fixed cut-off. In particular the quantities 

g (x) = <^\al j (x)a ij (x)\^ > (21) 

oo „i I n \ n 

= ^2 dxi... dx n 5 x m - 1 ^2 6(x m - x)\f n \ 2 , 

n=2^° \m=l / m=l 

< n > = [ dx g(x) , (22) 
Jo 

1 = dx xg(x) (momentum sum rule) , (23) 
Jo 

give the the number of gluons g(x) with momentum fraction between x and x + dx in a bound- 
state and the average number of gluons < n >. These quantities are helpful in trying to match 
the 1 + 1-dimensional eigentates forming represenations of discrete reduced symmetries with the 
representations J of the full spatial rotation group SO(2) in 2 + 1 dimensions. In addition they 
suggest appropriate analytic approximations. 



S 



3.2 Solutions. 



The numerical solutions for the reduced D = 3 theory are given in Table 1 and illustrated on Fig- 
ure 2 (some data for fixed K are tabulated in ref.Q). As in refs.|| |6| we note that a spectacularly- 
accurate valence gluon approximation is manifested in the light glueballs described by (|l~6|). The 
Coulomb potential, which includes the 2nd term in ( |i~6| ) as well as the -E-integral, dominates since 
it is the only positive definite and singular amplitude. Both the annihilation channel (D) and pair 
production (F) amplitudes take either sign with roughly equal probability, leading to much cance- 
lation. Physically they are suppressed because the intermediate A + 'particle' is non-propagating 
in light-cone time x + , so couples unfavourably to transverse gluon pairs in the boundstate. The 
same reasoning may be applied to the good valence approximation for QCD 2 with massive quarks 
at finite N, which has been studied in the same formalism [f[7f |. Therefore the latter is probably 
not a very good guide to the accuracy of valence approximations for quarks in higher dimensional 
QCD, where pair production can proceed via transverse gluons. 

In labeling the states with their (P\,C) quantum numbers, a technical problem arises in de- 
termining Pi, which has the effect x + «-> x~~ on vectors. This symmetry is broken when the 
momentum fractions x are discretised, since this is equivalent to (anti)periodic conditions on x~ . 
Under Pi, momentum fractions of n on-shell partons transform as |l7j 

x m — ► - ™ — , (24) 

Xm 2^m'=l Xm , 

which in general is not of the form (integer) /if; the exception is n = 2, in which case it has the 
effect x\ <-> X2- Also it is not obvious that eq(|24|) represents Pi in the interacting theory, although 
we found it to be at least an approximate symmetry of the states exhibited in this paper. |^| There 
are various quantitative and qualitative ways of assessing by inspection the effect of Pi as given 
by ( |24| ) on eigenfunctions at fixed cut-off K. The qualitative one we found particularly useful was 
to note that in general Pi transforms Fock states where one gluon carries most of the longitudinal 
momentum to Fock states where the momentum is shared evenly between all the gluons. 

A striking property of fig. 2 is the almost linear M versus < n > trajectory, which continues 
much further than we have shown and which is approximately repeated before soon becoming 
diluted by pair production effects at higher mass. Such a trajectory represents 'radial excitation' 
of the glueball flux loop in the sense that the mass of a state is increased by adding gluons and 
their attendant flux lines to the ring. The leading radial trajectory shows featureless structure 
functions (Figure 3) peaked at 1/ < n >, while the higher trajectory exhibits further oscillatory 
3 After submission of this paper, similar work appeared for N = 3 |E9], where a fuller discussion of parity can be found. 
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behaviour analogous to higher 'angular momentum' states. For the latter there is an increased 
probability of finding an asymmetrical sharing of the momentum between gluons in the glueball. 

It is interesting to note that g(x) remains large for x ~ in general, as a result of the massless 
limit m —* 0. In particular the approximately two-gluon states are rather like cosine wavef unctions, 
which was also found for the low- lying mesons in the 't Hooft model at small quark mass [17|. In 



the latter case the massless groundstate meson formed from massless quark and anti-quark had a 
constant wavefunction in phase space f2(x, 1 — x) = const.. In the present case the CP\ = ++ 
groundstate wavefunction is approximately gluonium 

|g >~ / dx f ^ X, t~ x) Tr[at(x)at(l - x))\0 > (25) 
Jo N 

The wavefunction f%(x, 1 — x) = g(x) is not quite constant due to the extra statistics-related 
contribution it j\^/x\X2 to its mass in ([jl]), suppressing it at the endpoints; this gives the glueball 
spectrum its mass gap. The gluonium approximation, where all f n are set to zero except for fa, is 
equivalent to the problem of scalar fundamental representation matter [18], with g 2 — > 2g 2 . The 
endpoint behaviour is ji ~ x a for x — > where 

acot(7r(- - a)) = ^— j , (26) 

which is important for determining the form of 1/K corrections (see appendix). Note that for 
m — > one has $2{x) — ► const > as x — > 0. By arguments analogous to those for quarks at finite 



N [17 1, we expect this endpoint behaviour to hold true beyond the gluonium approximation. 

On the leading radial trajectory the wavefunctions are well-described by constant ones in phase 
space for most momentum fractions x. Only near x ~ 0, corresponding to large x~ -separations, 
do interactions become appreciable. Therefore one can obtain an estimate for those masses by 
evaluating M 2 in a sub-basis of phase space wavefunctions (PSW) defined by 

ipn ■ fn(%l, ...x n )= const. , f m = for m ^ n . (27) 

A straightforward calculation using ( |l~6[ ) shows that M 2 is diagonal; 

M^n >= n(n "f g2i V > (28) 
g(x) = n{n - - x) n - 2 , (29) 



in fair agreement with the numerical data (table. 1). Only the 2nd term on the r.h.s. of eq. fllq) , 
7r l^y]x\Xii which is left over from partial cancellation of the processes involving radiation of zero- 
momentum A + (figs.l(a)(b)), contributes to the eigenvalue; it is the principal value of the E 
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integral. Since M ~ n asymptotically, one could regard it as the contribution to a constituent 
gluon mass due to longitudinal interactions, or more appropriately as the groundstate energy 
stored in the flux string between each pair of neighbouring gluons on the flux ring. Practically 
speaking, the PSW approximation works well because the lowest radial trajectory starts quite a 
bit lower than the complex higher trajectories. 

3.3 Comparison with QCD 3 . 

Recent lattice Monte Carlo data for SU(3) glueball masses extrapolated to the continuum limit in 
three dimensions § is given in table 1 for comparison. Finite lattice spacing errors are significantly 
reduced in three dimensions, so these results can be considered quite reliable. The classification 
is J PC i where parity P in two space dimensions is taken to be reflection about one axis (e.g. 

). Particles of non-zero angular momentum J should be degenerate in P- 
doublets | J > ±| — J >. Although we have no way of telling the true strength of breaking 
of SO (2) rotational symmetry by the dimensionally reduced theory, which reduces it to a Z2 
subgroup, we will attempt to match our states with J PC labels. Firstly we take P = P\. Under 
180° rotations x 1 — > — x 1 and x 2 — > — x 2 , but although cj) is a 1 + 1-dimensional scalar, it is the 
component of a 2 + 1-dimensional vector, so we will assume this induces P\ and (j) — > —<p. Since 
all the states we study are found to be invariant under 

O : <fiij -> <f>ji, (30) 

which reverses the ordering of gluons around a flux ring, this implies that CP\ = (— l)'" 7 ' mod2 ; 
which rules out half the possible states in three dimensions from the very beginning. In fact it 
removes precisely one of the states from each P-doublet. (We should point out however that 
states odd under O can of course be constructed, but that they lie consistently much higher in the 
spectrum.) One is still left with the problem of spin labelling within the set of even and odd spins 
separately, because in general there will be mixing within each set, i.e. spin mixes with spin 
2 etc.. We took the spin label for states of given (P\,C) on a radial trajectory to be the lowest 
which has not appeared with the same (Pi, C) on a lower energy trajectory. 

Having truncated to zero transverse momentum, the best one could hope for would be to 
reproduce the correct level ordering of the higher dimensional theory. The mass ratio scale gets 
expanded in the reduced theory, which one could ascribe to the level mixing as a result of SO(2)- 
breaking. It is tempting, in this gauge and in this large-iV light-front quantisation scheme, to 
describe the data in terms of a constituent gluon picture. However, in at least one case this is 
in clear conflict with the (N = 3) lattice data: a 0~~ state could be interpreted as an orbital 
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excitation of a 3-gluon groundstate l" 1 ; on the lattice the O is the lower state since it couples 
to the combination Up — Up of one Wilson plaquette Up, while l" 1 requires a longer lattice Wilson 
loop. In fact table 1 shows that the (N = oo) dimensionally reduced spectra in the C = —1 sector 
predict lighter spin 1 glueballs, in disagreement with the 3D lattice. The ordering in the C = +1 
sector is quite good by comparison. Spin l's being lower than spin 2's may be due to the mixing 
of even spins and odd spins amongst themselves since spin 1 is pushed down by mixing with spin 

3 etc. Due to the dominant valence content of wavefunctions, it is easy to assess the effects of a 
non-zero mass m for gluons. Although this would act to shrink the mass ratio scale and push up 
the spin l's, it could not produce a lower-lying O . 

Lattice results also exist for SU(2) when C = +1 ||, although states with C = — 1 do not 
appear since Wilson loops are unoriented. We also therefore tried to assess the error in extrap- 
olation to N = oo by fitting to the form = Mn — const. /N 2 for example. The variation is 
actually quite small and we do not think that finite iV corrections are the major source of error, 
relatively speaking, but that the discrepancies we have found are probably largely of a kinematic 
nature due to the severity of the truncation and hence the smallness of the residual Z<i symmetry. 
Despite some agreement of the data, our P and especially J assignments must be considered as 
heuristic since mixing is inevitable. In many ways the reduced model of QCD 4 , which is the case 
of physical interest, is less ambiguous. 

4 Reduced QCD 4 . 

4.1 Boundstate Equations. 

We now consider the reduced model for D = 4, for which there are two transverse polarizations. 
Since all particles are restricted to move in one space dimension it is especially convenient to work 
in the helicity basis (c.f.eq(|l3|)) 

4 = ^ ■ (31) 

Then 

t*nn 



h = J dk {a l +lj {k)a +ij (k) - aLyCAOa-tfC*)} ( 32 ) 

measures the helicity of a state. The glueball mass eigenstates will have definite helicity, since 
it is a symmetry of Sr, and fall into degenerate pairs of opposite h. In general of course, these 
will not exhibit the further degeneracies required to make up SO(3) representations, a penalty of 
throwing away all transverse momenta. 
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In the helicity basis the light-front momentum and energy (pi), (10) read 



f°° + + 

P+ = / dk fc{a] Hi (fc)a +ii (fc) + al^Ki^)} , (33) 



P~ = T {a ] +lJ (k)a +%J (k) + alyWa-yW} 

<7 2 f°° dkidk 2 dk 3 dk4 ,., nI , , . . + . . + . . . 

+ — / 7= , j , 0(fcl + K 2 - fc 3 - K4){(-D Kl)«2,K3 + («3ja_ {k^CL-U ( fci )a+/ fc (fe 2 ) 

o7T Jo \/k\k 2 k 3 k^ 

+( J D[fcl,A^,/c3] - #[fci,&2,&3] - 2t)a 1 fii (A;3)a^ fc ( fc 4)a+ l j(A;i)a-jfc(fc2) 

+ (i - E[ki,k 2 , k ? X)a\ ij (k^)a ] +jk (k i )a + ii(ki)a + i k {k 2 )} 

+8(kx + k 2 + k 3 - h){(F[ki,k 2 , h] - t)[a\ ij (ki)ai jk (k2)a\ kl (k 3 )a +i i(k 4 ,) 

+a+i j (hW_ jk (k 2 )a j +kl (k 1 )a +i i(k 4 ) + a 1 fi ;(/c 4 )a +ij (^i)a-jfe(^2)a+fcK fc 3) 

+a+ii(^4)a+ii(^3)a-ifc(^2)a+fc/(fci)] 

+ fe 2 , fc 3 ] + t)[ai ij (k 3 )al_ jk (k 2 )a\. ia {k 1 )a-ufa) 

+a+ij(h)a i _ jk (k 2 )a i _ kl (k 3 )a^ i i(k4) + a^A^a+ii (ksja+jki^a-kiih) 

+a\ii{k i )a^ ij {ki)a +jk {k 2 )a +k i(k 3 )]} 

+ (a+ <-> o_) . (34) 

The same set of elementary processes occur (fig.l), with helicity conserved at each vertex. A set 
of integral equations similar to ( |i~6| ) may be derived; in the following the polarization index is now 
taken as p = ±, with the properties 5-\ = <5 + + = 1, 5 ++ = 0, etc.. 

■^fftn-Pn = f P v--p n (xi,x 2 ,...,x n ) + ^=f Pl ... Pn (x l ,x 2 ,...,x n ) 

+ / —, j {E[x 1 ,x 2 ,y]f pi ... pn (x 1 ,x 2 ,,...,x n ) 

Jo ^x\x 2 y[x\ + x 2 - y) 

+ (D[xi,x 2 ,y] +t)S PlP2 f p pp3 ... pn (y,x 1 + x 2 - y, . . . , x n ) 
+ ( 4t6P p 2 ~ 3t ~ E [ x i^ x 2,y])f P i- Pn (y^ x i + x 2 -y,.. .,x n )} 
+ / dy 



Jo Ay/xiyz(xi - y- z) 



{(F[x!,-y,-z] -t sgn(pi)sgn(p))/ p / pp2 ... pn (y,z,xi - y - z,x 2 , . . . ,x n ) 
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+ ( F [-z, -V, x i\ + t s S^(Pi)sgn(p))f p ppi ... pn (y,z,xi-y- z,x 2 ,..., x n )} 

+ . , , l , ; A (.F[x 3 , x 2 , xi] + 1 sgn(pi)sgn(p 2 )) x 

4^/xiX2X 3 (3;i + x 2 + x 3 ) 

fpiP4-Pn ( X l + X2 + X 3 , X 4 , • • • , S„) 

+ (F[x 1 ,x 2 ,x 3 ] +t sgn(p 2 )sgn(p 3 ))dp 1 p 2 f P3 ... Pn (x 1 + x 2 + x 3 ,x 4 ,. . . ,a%)} 

+ cyclic permutations of (x\,x 2 , . . . , x n ) . (35) 

The mass renormalisation proceeds in exactly the same fashion — the only new ingredient is the 
normal ordering of the t-terms which gives a further logarithmically divergent contribution to 
mf nd — and it is in the form ( poj ) that we discretise the momentum fractions x for numerical 
computation. 

The variable t is considered as a free parameter, which we will adjust to obtain a physically 
reasonable spectrum. In particular we find that it must be positive in order for the groundstate 
glueball to be helicity zero^J Unlike in QCD 3 , gluon pair production can now proceed directly from 
a contact term fig. 1(c). While it is therefore less suppressed, for a range of t around zero there 
is still a discernable valence gluon content for light glueballs. If t is too large however, tachyons 
appear and eventually there appears to be a phase transition; the phase structure is discussed in 
more detail in section 5. As in section 3, we primarily study the massless case m = 0. However 
we will find once again that there is an expansion of the mass ratio scale in the reduced theory 
relative to that obtained from the higher dimensional lattice results, which can be mitigated by 
including a constituent mass m by hand. 

From the wavefunctions (^) we also calculate both unpolarized 

n( , <V\al ij (x)a +ij (x) + al ij (x)a- ij (x) \V > 

G(x) = J - — - — J - (36) 

< * * > 



and polarized 



AG{X) = ^ ( 3? ) 



4r 



This result was also found for the similar transverse plaquette coupling in the light-front hamiltonian lattice gauge 
theory of ref.[|). 
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structure functions. As a consequence of no transverse momenta, there is a trivially satisfied 
helicity sum rule 

/ dx AG(x) = h . (38) 
J o 

In higher dimensions these structure functions would also be integrated over transverse momenta 
up to some scale Q, the evolution with large Q being calculable within perturbative QCD once the 
function is specified at a given scale. In the present case Q = effectively, and we are performing 
a non-perturbative calculation of (some of) the long-distance properties of boundstates in QCD 4 
scarcely known for any hadron, least of all glueballs. 

4.2 Solutions. 

Results for masses of the seven lightest glueball states at t = 0.15 (m = 0) are given in Table 2 
and illustrated on Figure 4 classified by l^l^ 1 . For massless gluons m = we are forced to choose 
a relatively small t in order to avoid tachyons, so it is the longitudinal interactions which are still 
dominant. As expected, the valence approximation is less accurate than for the D = 3 reduced 
model, but is nevertheless manifest. All the states shown are groundstates with respect to their 
valence gluon content and, with the exception of , exhibit featureless unpolarized structure 
functions (Figure 5) peaked at 1/ < n >. The x —* behaviour is still governed by eq.(|26|). The 
negative Pi-parity state , like the for D = 3, has strong longitudinal interactions, which 
is why it appears relatively heavy. We attribute the persistence of valence gluon dominance to 
the energy required to add extra flux lines (one per gluon) to the ring, since the study of purely 
longitudinal interactions in section 3 showed that each contributes roughly ~ g^NTr/8 to the 
mass. 

Once more, it is possible to estimate the light spectrum in the PSW approximation by setting 

f Pl ... Pn (xi, ... ,x n ) = const. (39) 

for a given set of gluon helicities p±,...,p n . This provides a sub-basis of wavefunctions tp pi ... Pn for 
the 'radial' excitations of the groundstate in each \h\ sector (Pi = +). The integrals appearing 
in the Bethe-Salpeter equations (|35| ) can be done exactly for ( |39"[) resulting in a finite M 2 matrix 
for finite maximum n. Note that M 2 is now no longer diagonal in each sector due to the t-terms. 
Diagonalising in the sub-basis of PSW's up to n = 6 gives the estimates shown as Mp SW in table 2. 
As an example we show the M 2 sub-matrix in the PSW basis {tp^ , V'+h , ip-\ I > > VM-+h } : 
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f ^(l-i) 2vrV6t -4irV3t \ 

2vr^/6t 3vr 2 (l-t/2) 

-£iry/?,t 2*Jt 3vr 2 (l-2f) 

+ 

\ 3vr 2 / 



(40) 



Polarized structure functions of /i = glueballs are identically zero while that of the \h\ = 2 
state follows closely the unpolarized function because the non-valence content is small. The spin 
structure of the \h\ = 1 state is perhaps the most interesting since it is the nearest thing to a 
'nucleon' in the spectrum. Its larger non-valence content means that the unpolarized function is 
peaked at noticeably smaller x than the polarized one, gluon pairs, which are preferably created at 
small x, having to be produced with opposite helicity. In Figure 6 we plot the helicity asymmetry 
A = AG/G for the l" 1 glueball. The PSW approximation would give A(x) = 1/3. At small x one 
finds A(x) < 1/3 so the helicity tends to be more disordered, while at large x one has A(x) > 1/3. 



In real nucleons one expects A —* 1 as x — > 1 [19| (helicity retention) i.e. a parton carrying most 



of the momentum also carries most of the helicity. Although this question has presumably not 
been asked of glueballs before, it does not seem to happen in the reduced model for t = 0.15, at 
cut-off up to K = 13 at least, although the x = 1 intercept of A does increase with t. Perhaps 
this an indication that the strength of transverse interactions we are forced to consider at m = 
are too small to accurately reflect the higher dimensional theory. 

4.3 Comparison with QCD 4 . 

Lattice Monte Carlo simulations of the (quenched) glueball spectrum have been relatively rare and 
are still crude for all but the lightest state (see e.g. refs.|2C, ^]). We used the SU(3) results of 



ref.[22| for comparison with the large- N reduced model since they include a range of light states. 
However that computation was performed at fixed lattice spacing (j3 = 6.4 on a 32 4 lattice) so 
one must bear in mind that significant variations in mass ratios are possible when extrapolated 
to the four-dimensional continuum limit; furthermore, no 'radial' excitations were computed. It is 
generally believed that N = 3 is large in the sense that there is little qualitative difference between 
SU(2) and SU(3) spectra. Lorentz invariance was restored in the lattice results to the extent that 
a full J PC classification was possible, but this is clearly not the case in the reduced model, where 
are the exact quantum numbers. As in the D = 3 reduced model, some of the components 
of given light spin J are pushed higher in spectrum, although each light spin J does seem to be 
represented by its largest \J Z \ values (compare helicity h with spin J in Table 2). In particular the 
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reduced model does not seem to cope well with states odd under O : <pij — > <fiji, which are needed 
to make up components with C = — (— l)l Jz l mod2 (e.g. J z = ±1 for J PC = 2 ++ ). The results of 
Table 2 also show that mixing between given seems to prevent the degeneracy with h = 2 of 
an h = component, which mixes with excitations of the groundstate scalar glueball ++ . Finally 
we note that Pi alone seems to be closer to three dimensional parity P, rather than Pi combined 
with (j) p — > —cj)p, since the parity of the lightest spin 1 would be incorrect in the latter case. 

Assuming that the h assignments of the reduced model are a good guide to the J assignments in 
four dimensions, the level ordering is in qualitative agreement with the lattice result, although once 
again the mass ratio scale is expanded. Our model (with m = 0) prefers a heavy '", due to strong 
longitudinal interactions, which we have treated exactly, while the four-dimensional lattice would 
prefer it to be lighter. It is quite possible that the full effects of transverse interactions, of which we 
have only retained a remnant, and transverse kinematics shrink the mass ratio scale and make the 
h lighter. To illustrate this possibilty one can include a constituent mass m by hand, which one 
might regard as generated by the transverse dynamics that the reduced model neglects. In Table 
2 we have also given the spectrum at K = 12 for a set of parameters m 2 = 2g 2 N/ir,t = 0.5. The 
non-zero gluon mass allows us to choose a value for t comparable to the longitudinal interaction 
strength. The comparison with the (N = 3) four-dimensional lattice data shows much closer 
agreement, although no amount of fudging can cause the degeneracies necessary for restoration of 
SO(3) rotational symmetries; the transverse kinematics are clearly needed for this to happen. The 
main effect on the structure functions is to suppress small x partons (c.f. eq.^)), the polarization 
asymmetry A also being modified slightly (fig. 6). 

5 Phase Structure and Matrix Models. 

In principle there are two adjustable parameters in the reduced model for QCD 4 , the mass m and 
quartic self-coupling t. In this section, by deriving bounds on M 2 and quantifying independence, 
we suggest a phase diagram in the (m 2 ,t) plane and note the analogy with matrix models of 
non-critical string theory. As is well known, the planar Feynman diagrams for the field theory of 
a hermitian matrix cpij with action 



S = / d d xTr 



^ AN 



(41) 



define a certain type of random surface embedded in D dimensions [23], essentially a discretisation 



of Polyakov's continuum path integral formulation [24]. Such theories have been analysed for d = 2 



and N = co va light-front formalism [13, 14 1 in very much the same way as the two-dimensioanl 
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adjoint gauge theories. The action ( |4l| ) is merely ungauged. The critical value of X/fi beyond 
which perturbation theory diverges — individual normal-ordered planar diagrams are UV finite 
- manifests itself in light-front formalism (without zero modes) as the divergence M 2 — ► — oo in 
the spectrum. Such a transition^] occurs for the reduced model of QCD 4 also, separating a parton 
phase with < n > finite from a string phase with < n >= oo. 

Firstly, one can prove that the spectrum is unbounded below if m 2 < by using the variational 
wavefunction 

/ dx Tr[a\(x)al(l - x)]\0 > . (42) 
Jo 

Secondly, that the spectrum is unbounded from below if t is sufficiently large can be proven with 
the (discretised) variational state 

TT[al(l/K)al(l/K)al{l/K)al{l/K)...al(l/K)al(l/K)}\0 > . (43) 

Therefore the phase diagram is probably as indicated in Figure 7. In principle it is possible that 
the point A coincides with m = t = and that the spectrum is unbounded below for all t > 
when m = 0. The l/^f -extrapolation curve of the groundstate M 2 shown in the appendix certainly 
seems to be dropping rapidly. However, although we could not prove it, our experience suggests 
that unboundedness is unlikely since we could not find a variational wavefunction which proved 
that M 2 = — oo for all t > at m = — it is certain that M 2 is positive at t = m = — while 



the matrix model analysis makes definite predictions for the K dependence [14], which can be 
manifested at quite small K, for the divergences of the string phase; 

< n >~ K , M 2 K 2 . (44) 

Our measurements at t = 0.15 showed no sign of a linear divergence of < n >, rather it seemed 
to converge to ~ 2.06. However, at t = 0.45 for example we find the -ff-dependence of eq. (|44]) , 
suggesting that this point is already in the string phase. The precise details at the transitional 
values of t are difficult to determine at the small values of K with which we are forced to work here, 
but it would be interesting to know whether there exists a non-critical string regime analogous 



to that identified in refs.[13, 14]. Another question that we are unable to answer at this stage 
concerns the role of k + = modes on the phase structure, since they probably have to be taken into 
consideration in the string phase for a faithful treatment of the original field theory (g). Finally 
we note that if our conclusion concerning boundedness were wrong for m = 0, one would be forced 
to introduce a sufficiently large gluon mass m for given t in order to stabilize the spectrum. 



The possibility of phase transition in the adjoint gauge theory when contact self-interaction terms are added was 
already suggested in refe.[p5[ |). 
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6 Conclusions. 



Two-dimensional non-abelian gauge theories have proved to be usefully tractable, especially in 
light-front formalism and the large- N limit. In this paper we have attempted to incorporate some 
of the transverse degrees of freedom of pure gauge theory in a two-dimensional context. In this way 
one can investigate the boundstate properties of physical gluons and spin while postponing many 
of the subtle problems of renormalisation which arise in higher dimensional light-front formalism. 
Although our asymmetrical truncation of phase space may seem an arbitrary approximation, the 
highest penalty paid appears to be the lack of Lorentz multiplet degeneracies in the glueball spec- 
trum. Comparison with the lattice Monte Carlo data for both three and four dimensions showed 
nevertheless that the dimensionally reduced model always produces some components of each light 
Lorentz multiplet, and in approximately the same level order. The rogue components which do 
not appear in the light spectrum are states anti-symmetric under the symmetry transformation 
Aij — > Aji of the gauge potential; this fact deserves a better understanding. Although one cannot 
be sure of the significance of the mechanism in higher dimensions, the two-dimensional model 
simply generates from the longitudinal Coulomb potential what amounts to a constituent gluon 
mass structure for the spectrum. This is due to the non-zero groundstate energy for the flux 
string between neighboring gluons in the flux loop. It is weaker than a true gluon mass and so 
cannot prevent the wavefunction from spreading in x-space when the gluon mass m = 0, although 
there is a mass gap in the spectrum. It was also the origin of a good valence approximation for 
the D = 4 reduced model. Such an approximation is even more evident for the D = 3 model 
in which case pair production must proceed via a non-propagating longitudinally polarized A + 
field. An improvement in the comparison with the higher dimensional spectrum was noticed when 
an explicit gluon mass m is included by hand — this might be neccessary anyway to bound the 
spectrum from below — perhaps compensating for neglected transverse dynamics. The light-front 
structure functions exhibited physically reasonable behaviour in our chosen region of coupling 
constant space, being quite close to those of a phase space approximation. More interesting will 
be their content at finite N, which will govern the decay and formation of glueballs, and when 
quarks are included. 

Although many qualitative features of higher dimensional QCD are reproduced by the reduced 
models, we are sceptical of its use as the basis for a systematically improvable approximation 
to Lorentz-invariant quantities at the quantitative level as transverse momenta are added. The 
limit of zero transverse momentum is a very singular one in the quantum theory |15[ |, which we 
have avoided by reducing at the classical level and treating subsequent undetermined couplings 
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phenomenologically. Of course we would be happy to be proved wrong in this scepticism! Rather 
we consider the transverse lattice approach advocated in ref.|J to be the most promising for 
quantitative work — a study is currently in progress — whose mathematical structure is virtually 
isomorphic to the field theories studied in this paper, although the adjoint fields have a different 
interpretation. Nearly co-linear phenomena in which transverse momenta are small compared to a 
typical scale, such as high-energy diffractive scattering, may hold the best opportunity for applying 
the reduced models quantitatively to a bona fide four-dimensional problem. With this in mind, the 
scattering amplitudes involving the reduced glueballs and mesons [3C] are presently under study. 
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APPENDIX. 



The Lanczos Algorithm. 

The Lanczos algorithm [p6| , p7| is an iteration method which may be used to determine numer- 
ically the eigenvalues and corresponding eigenstates of a finite linear operator. We may summarise 
the essential ideas as follows: Given a symmetric eigenvalue problem Au = Xu, we begin by choos- 
ing a normalised vector u 1; and setting the constant 61 to zero. A sequence of constants a n ,b n 
and orthonormal vectors u n are then generated by the following recursion relations; 

v n+ i = Au n - 6 n tt n _i, 
a n = v n+ i ■ u n , 

v n+l = v n+l ~ Q>n u m 



u n +i = v' n+1 /b n+1 . 

It is easy to show that the matrix representation of the operator A with respect to the orthonormal 
vectors {u n } is tridiagonal, with a n and b n being the diagonal and codiagonal entries respectively. 
After a given number of iterations, the eigenvalues and eigenvectors of this tridiagonal matrix 
approximate those of the original operator, with a progressive improvement in precision as the 
iteration proceeds. A particular feature of this method is that the extreme eigenvalues are readily 
obtained to high precision after a comparatively small number of iterations. This makes the 
algorithm especially useful if one is only interested in studying the low-lying states. 

In our investigations, the finite linear operators which need to be diagonalised arise from the 
discretised versions of the integral equations (|l6|) and (|35| ) , in which only odd[] momentum fractions 

Xi £ {1/K,3/K,5/K, . . .} are allowed, coupled with the constraint x\ -) Yx n = 1. This enables 

us to replace the functions f pl ... Pn (xi, . . . , x n ) with a finite number of components 

9pi...p„[mi, • • • , m-ii) — ^j -^n-i Jpi—Pn\ 1 ' ' ' 1 

where the odd integers rrii satisfy m\ + • • • + m n = K. Note that this cutoff regulates both 
momentum and the number of particles. Thus the linear action of the integral equation (35) on 



the space of functions f Pl ... Pn may be represented in the discrete case as a finite constant matrix 
whose size is determined by the number of independent components g pi ___ Pn (mi, . . . ,m n ). For 
example, if we set K = 4, and consider the zero helicity sector of the reduced 3 + 1 theory, 



Choosing odd multiples gives rise to a faster rate of convergence towards the continuum limit since the endpoint 
regions of the x-interval (0, 1) are better sampled. 
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the only relevant components are (l,3),g_| (3, 1),<?+h (1) 1; 1> 1)> an d Sh I (1> 1) 1> !)■ The 

corresponding 4x4 symmetric matrix representing the discretised version of the linear action (|35|) 
has the explicit form 



/ 2tl , 16 _ 4y. , 16 _16_4 / _4_y. 8_y. \ 

/ 3 3 + 3 3 3 yr \ 

_16 _ 4 + 2tt i 16 _ 4 + , 16 _4_ f 8 f 

3 3 6 ^3 ~ t ~ 3 3 ^ 3 V6 

v¥ v¥ 4vr + 16*-4t ^ 



(45) 



V -76* 75* 4vr + 16x-16ty 

where cc = -Km 2 / g 2 N . Diagonalising the above matrix gives the eigenvalues M 2 in units g 2 N/ir. 

One strategy which may be used to estimate the continuum (i.e. K — > oo) limit of the theory 
is to calculate these matrices explicitly for increasing values of K, and then to extrapolate the 
resulting sequence of eigenvalues. Unfortunately, this procedure is severely limited by the rapid 
growth in the number of independent components g Pl ... p „(mi, . . . ,m n ) as K is steadily increased, 
making it sometimes difficult to provide reliable estimates of the continuum behaviour of the 
theory. Approximation methods are then necessary if further progress is to be made. 

As in a previous investigation pip , we used Sun and Dec-alpha workstations together with 
Mathematica to implement the Lanczos algorithm outlined earlier — this way is slow but easy 
to do. This technique enabled the authors to investigate the cases K < 13 in the reduced 3 + 1 
theory and K < 22 for 2 + 1. Typically, the algorithm requires a judicious choice of initial vector 
to ensure that the convergence of the low- lying eigenvalues is well behaved and rapid. Owing to 
the valence-like structure of the low-lying states, a perfectly adequate initial state consists of two 
or three valence gluons (depending on whether the helicity is even or odd respectively). It is also 
important that it have definite symmetry under the operators which commute with the light-cone 
hamiltonian since, as the iteration proceeds, each of the new vectors formed by the action of the 
light-cone hamiltonian will possess the same quantum numbers as the initial vector. This helps 
restrict the size of the vector space one needs to work with. We remark that for K = 12, the ++ 
sector consists of approximately 900 linearly independent vectors, while for K = 13, the size of 
the l" 1 sector is a little over 1800. 

In the context of discretised light-cone quantisation, exploiting the advantages offered by the 
Lanczos algorithm depends crucially on how one implements the method. Typically, a large 
number of vectors are produced at each stage in the iteration, and so, as the iteration proceeds, 
the light-cone hamiltonian will repeatedly act on the same vectors. However, this action is far 
from trivial, and so the procedure rapidly runs into difficulties due to the enormous growth in 
the number of vectors that need to be dealt with. Therefore, it was necessary to introduce a 
system of 'pointers' which connect a set of existing vectors with its image under the action of the 
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hamiltonian. This conveniently eliminates the need to perform the same calculations repeatedly, 

and dramatically improves the viability of the method. 

Extrapolation. 

A fool-proof method of extrapolating the results of discretised light-cone quantisation to the 
continuum limit is still lacking. The extrapolation of finite- K calculations to the continuum limit 
1/K = has been discussed in refs.[[l7]]. Errors due to discretisation of the integrals in equations 
such as ( |i~6| ) begin at 0(1/ K) due to the Coulomb pole. In addition, endpoint errors lead to 
dependence upon K~ a (c.f. eq(p6|)). Thus for the case m = one expects the corrections to form 
a series in 1/K and one can fit n data points taken at different K to a polynomial (Richardson 
extrapolation) 

with ao the continuum limit. a n ^i/ K n ~ l evaluated at the lowest value of K employed gives 
an order of magnitude estimate for the error in extrapolating . These are the errors shown 
in the tables, which are probably underestimated. We illustrate this procedure below with the 
masses from the D = 4 reduced model (section 4) at t = 0.15, since the extrapolation is far from 
straightforward especially for the groundstate in this case. Given the suppression of non-valence 
gluons, one can increase the maximum value of K attainable numerically by truncating the Fock 
space to a few gluons (Tamm-Dancoff approximation), thus enabling a more reliable extrapolation 
to K = do. Figure 8 illustrates the polynomial fits for K < 12, 13 for all states shown in table 
2, and also includes the fit for the groundstate obtained from data for K < 18 by using the 
2— , 4— , 6-gluon sector only (chain line). 
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FIGURE and TABLE CAPTIONS 



Table 1 - M is the K = oo extrapolated mass in units of y/Ng 2 /n based upon exact diagonal- 
isations up to K = 22. < n > is calculated for illustration from the theory truncated to the 
2-,4-,6-gluon sector only for K = 20, or to the 3-, 5- gluon sector for K = 19. The heuristic 
J assignments are discussed in text. The lattice SU(S) masses Mgu^ have groundstate ++ 
normalised to M and a la statistical error estimate (ours). 

Table 2 - Exact M 2 in units of Ng 2 /ir for various K at m = and t = 0.15; the continuum limit; 
< n > at K = 12; = M 2 at K = 12 for m 2 = 2g 2 N/ir and i = 0.5; < n > for the latter. 
Also given are mass squared ratios M 2 to the groundstate. Lattice SU (3) mass squared ratios to 
the groundstate M 2 ,^^ are taken from the calculation of ref.|2^] at fixed lattice spacing and are 
rounded to two figures here. 

Figure 1 - The elementary processes contributing at order g 2 in Sr: (a) self-energy; (b) longitu- 
dinal interactions; (c) transverse interactions. 

Figure 2 - Mass spectrum, in units of \/Ng 2 Jtt, extrapolated to the continuum limit K = oo 
and classified by P\,C. 

Figure 3 - Gluon structure functions g(x) of glueball mass eigenfunctions. Full lines refer to 
states on the lower 'radial' trajectory in fig. 2 and chain lines to the states on the upper trajectory, 
(a) < n >~ 2 ; (b) < n >~ 3 ;(c) < n >~ 4 ; (d) < n >~ 5. 

Figure 4 - Mass squared spectrum in units of Ng 2 /ir at m = and t = 0.15, extrapolated to the 
continuum limit K = oo and classified by l/?,^ 1 '^ (non-zero \h\ are doublets). 

Figure 5 - Gluon structure functions G{x) for m = and t = 0.15 

Figure 6 - Polarization Asymmetry A(x) = A.G/G of the l" 1 glueball at K = 13. Dotted line 
for m = and t = 0.15, full line for m 2 = 2g 2 N/-K and t = 0.5. 

Figure 7 - Suggested (m 2 ,t) phase diagram. The shaded region has unbounded M 2 — > — oo, 
while the theory is non-tachyonic in a neighbourhood of the postive m 2 -axis. 

Figure 8 - Polynomial fits to M 2 data (units of g 2 N/-K) at t = 0.15 and m=0. 
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Table. 1 



REDUCED 




PiC(J) 


+ + (0) 


+ -(1) 


+ + (0) 


+ -(1) 


+ + (2) 


--(0) 


+ + (0) 


M 


2.1(1) 


3.4(1) 


4.7(1) 


5.9(1) 


6.3(1) 


6.6(1) 


7.1(1) 


< n > 


2.002 


3.002 


4.000 


4.990 


2.175 


3.085 


5.902 


M PSW 


2.2 


3.8 


5.4 


7.0 






8.6 


LATTICE 




J PC 


0++ 


0++ 


0~ 


2 ±+ 


0~ 


0++ 




M SU(3) 


2.1(1) 


3.2(1) 


3.2(1) 


3.5(3) 


4.0(2) 


4.0(2) 





Table.2 



REDUCED 




\h\* c 


0++ 


2++ 


1+- 


0++ 


0++ 


2++ 


o-+ 


valence 


a+a_ 


a±a± 


a±a±a,zf 




a+a+a—d- 


a±a±a±azf 


a_|_a_ 


K = 6 


3.12186 


4.28929 


8.78276 1 ) 


11.453 


14.502 


14.814 


16.1509 


K = 8 


3.24943 


4.43805 


9.18900 2 ) 


12.303 


15.813 


16.1555 


17.2769 


K = 10 


3.25875 


4.53335 


9.46014 3 ) 


12.835 


16.686 


17.0475 


18.0253 


K = 12 


3.25605 


4.59893 


9.65198 4 ) 


13.193 


17.309 


17.6828 


18.5584 


if = oo 


3.1(3) 


4.9(1) 


10.7(2) 


14.9(4) 


21.0(4) 


21.3(9) 


21.7(3) 


< n > 


2.062 


2.002 


3.031 4 ) 


4.011 


4.021 


4.025 


2.134 


Mp SW 


3.1 


5.7 


13.3 


17.0 


27.9 


28.6 




M 2 

±vi m 


14.56 


18.21 


37.29 4 ) 








35.98 


<n> 


2.042 


2.005 


3.027 4 ) 








2.022 


M 2 


1 


1.25 


2.56 








2.47 


LATTICE 




J PC 


0++ 


2++ 


1+- 








o-+ 




1 


1.5(1) 


1.9(2) 








1.5(2) 



VK = 7, = 9,^K = U^K = 13 
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Figure 1 
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